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ABSTRACT 

We  derive  a  formula  for  the  viscous  decay  of  long  Internal  solitary 
waves  that  propagate  Into  a  quiescent  fluid  In  a  two-layer  model.  The 
result  Is  analogous  to  Keulegan's  (1948)  formula  for  the  viscous  decay  of 
long  surface  waves.  The  requirement  that  the  fluid  ahead  of  the  wave  be 
quiescent  Is  Important,  and  we  show  experimentally  that  the  accuracy  of  the 
formula  decreases  significantly  If  the  Internal  waves  are  preceded  by 
faster-traveling  surface  waves. ^ 


1  Introduction  and  Analysis 


•> 


Solitary  water  waves  measured  In  laboratory  wave  tanks  are  known  to 
exhibit  appreciable  decay  due  to  viscous  boundary  layers  on  the  walls  of  the 
tank.  Keulegan  (1948)  developed  an.  approximate  theory  to  predict  this 
viscous  decay  of  surface  solitary  waves. 

Long  Internal  waves  resemble  long  surface  waves  In  may  respects,  and 
Figure  2  of  the  preceding  article  (hereafter  called  I)  demonstrates  that 
long  Internal  waves  also  experience  significant  viscous  dissipation.  If  the 
pycnocline  Is  thin,  then  the  dissipation  of  long  Internal  waves  Is  even 
stronger  than  that  of  long  surface  waves,  because  there  Is  an  additional 
viscous  boundary  layer  at  the  interface  of  the  two  fluids.  Thus,  to  account 
for  the  slow  viscous  dissipation  of  an  Internal  KdV  sollton.  It  is  necessary 
not  only  to  make  relatively  minor  changes  In  Keulegan's  analysis  of  the 
boundary  layers  on  the  walls  of  the  wave  tank,  but  also  to  account  for  the 
dissipation  at  the  Interface. 

The  purpose  of  this  paper  Is  to  study  the  viscous  dissipation  of  an 
Internal  KdV  sollton  as  It  propagates  down  a  laboratory  wave  tank.  The 
analysis  Is  based  on  the  original  work  of  Keulegan  (1948).  The  major 
assumptions  are  the  following. 

I)  There  are  two  fluids  of  different  densities  In  a  stably  stratified 
configuration  bounded  below  by  a  rigid  horizontal  bed  and  above  by  a 
free  surface  (see  Figure  1  of  I). 

II)  The  Inviscid  flow  Is  that  of  a  two-dimensional  Internal  wave  In  the 
form  of  a  KdV  sollton,  traveling  with  a  constant  speed  c. 
Consequently  we  may  write  a/at  ■  -c  a/  a^,  where  g  •  x-ct.  Moreover, 
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the  viscous  time- scale  Is  assumed  to  be  slow  In  comparison  with  the 
KdV  time-scale. 

Ill)  Nonlinear  effects  may  be  neglected  in  the  thin  viscous  boundary  layer, 
even  though  they  are  Important  In  establishing  the  Inviscid  flow. 

1v)  All  streamlines  may  be  cons1der«d  flat  and  horizontal  for  the  purpose 
of  the  viscous  analysis. 

v)  Only  shear  parallel  to  the  Interface  (which  Is  now  considered  flat)  Is 
of  consequence. 

v1)  Diffusion  of  species  Is  neglected,  as  Is  compressibility. 

With  these  assumptions,  one  may  show  that  the  larger  and  smaller  roots 

of 


(c2/g)2  -  (c2/g)(hi  +  h2)  +  Ahih2  ■  0  (1) 

define  the  propagation  speeds  of  the  long  surface  and  Internal  waves, 
respectively.  In  the  linear  Inviscid  limit  (see  I  or  Lamb,  1932,  section 
231).  Here  g  Is  the  constant  gravitational  force,  ^  ■  1-  Pj/p2 
dimensionless  density  difference,  and  hi,h2  are  the  upper  and  lower  fluid 
depths,  respectively.  We  assume  that  the  two  sets  of  roots  are  distinct, 
which  Is  guaranteed  If  either  a  ^  1  or  hi  ^  h2.  Let  ci  denote  a  solution  of 
(1)  corresponding  to  an  Internal  wave. 

The  linear  theory  also  defines  the  horizontal  velocities  In  the  two 
layers  In  terms  of  the  displacement  of  the  Interface.  In  dimensional 
variables,  the  velocities  of  the  lower  and  upper  fluids,  respectively,  for  a 
long  internal  wave  are 
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Mhere  n  Is  the  displacement  of  the  interface.  These  provide  the  "outer" 
inviscid  solution,  which  drives  the  viscous  boundary  layer  at  the  interface. 

The  KdV  equation  governs  the  evolution  on  the  next  time-scale  of  the 
long  Internal  wave  traveling  to  the  right.  The  dimensional  equation  Is 
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where 
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In  the  Boussinesq  limit  (a  ♦  0,  gA  finite),  ♦  gAh^h2/(hj  +  hj),  o  ♦  1, 
6  ♦>1  [h^^-  h"J],  -r  ♦  hjh^/e,  and  (3)  reduces  to  equation  (12)  of  I. 

As  a  further  check  on  the  coefficients  in  (3),  we  note  three  limits  in 
which  (1)  and  (3)  reduce  to  the  dimensional  form  of  the  KdV  equation  for 
surface  waves  on  a  homogeneous  fluid: 


(5h)-l/2li  *  li  4  liii  *  ^>4 

'  3t  8X  2  h  3x  6 


These  limits  are: 


1)  A 


~  ♦  (hj  ♦  h2) 


JL  m  - 5 -  i 

h2  h^  +  h2 


11)  Pi  ♦  0  ,  ^ 


111)  hi  -  0  . 
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Djordjevic  and  Redekopp  (1978)  also  derived  a  KdV  equation  for  long  Internal 
waves  without  the  Boussinesq  limit.  Their  results  differ  from  those  In  (3) 
except  In  the  limit  a  0,  because  they  assumed  a  rigid  lid  on  top,  whereas 
(3)  Is  based  on  having  a  free  surface  there.  The  boundary  condition  at  the 
upper  surface  is  unimportant  only  If  a.«  1. 

We  turn  now  to  the  problem  of  the  viscous  boundary  layer  at  the 
Interface.  In  dimensional  form,  the  viscous  equations  above  and  below  the 
dividing  streamline  are 


Ul) 


Ul) 


z  >  0 


(6) 

-  Ci  —  (U2  -  U2)  ■  V2  {U2  -  U2)  ,  2  <  0  , 

where  5  ■  x  -  Cit,  uj(5),  U2(c)  are  the  Inviscld  horizontal  velocities  In 
the  upper  and  lower  fluid,  u^(4,z),  U2(t,z)  are  the  solutions  of  the  viscous 
equations  In  the  two  layers,  and  vi,  v2  are  the  two  kinematic  viscosities. 
The  boundary  conditions  are  that 

Uj  ♦  Uj  as  z  ♦  ♦  •  , 

(7a) 


U2  U2  as  z  * 


and  at  z  *  0, 


Ui  •  U2  . 


(7b) 


where  wi,  u2  *re  the  two  dynamic  viscosities.  The  solutions  of  (6)  may  be 
written  in  the  form  of  the  Duhamel  integral  solution  of  the  heat  equation: 


Ui  -  Ui 


-  FiU  ♦  (c22/4vi  x^)]  exp(-  x^)  dx 


(8) 


U2  -  U2 
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F,C5  +  {czV4v2  x^)]  exp(-x^)  dx 


where 


Fi(0 


CUi(t)  -  U2(t)3  . 


(9) 


F2{0 


— ■ — — / — 7171  Cyi(4)  -  U2(5)3  » 
mi(v2)^'^  ♦  vzM 


and  we  require  that  u^(C).  U2(()  be  localized  enough  that  the  Integrals  In 
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(8)  converge. 


Next  we' need  an  expression  for  the  rate  of  energy  loss  due  to  viscous 
effects  at  the  Interface.  Consistent  with  assumption  (v)  above,  the  local 
dissipation  function  near  the  Interface  In  each  layer  Is 


j  “  1.2 


(10) 


The  total  rate  of  dissipation  due  to  the  boundary  layer  at  the  Interface  Is 
obtained  by  Integrating  ej  over  the  volume  of  that  boundary  layer. 
Integrating  once  by  parts  and  using  (6).  (7)  shows  that  the  rate  of  energy 
loss  due  to  the  boundary  layer  at  the  Interface  Is 


2  -0  ’  j  •  1  or  2  ,  (11) 

where  W  Is  the  total  width  of  the  channel.  Keulegan  (1948)  pointed  out  that 
(11)  represents  both  energy  dissipation  1n  the  boundary  layer  and  kinetic 
energy  left  behind  In  the  channel.  We  will  discuss  the  latter  effect  again 
in  Section  2. 

Energy  Is  also  lost  In  the  boundary  layers  on  the  walls  of  the  wave 
tank.  Keulegan  (1948)  showed  that  the  rate  of  energy  loss  In  each  of  these 
boundary  layers  Is 


■—  »  W  (ui -U2)[vj(8Uj/32)] 


dt 


Uo[)i(9u/3y)]  d(  , 


(12) 
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Mhere  az  Is  the  width  of  the  layer,  uo  Is  the  Inviscid  velocity  outside  the 
layer,  and  the  shear  stress  Is  evaluated  at  the  wall.  The  total  rate  of 
energy  loss  Is  obtained  by  combining  these  separate  contributions.  Including 
(11). 

We  now  refer  specifically  to  .a.soliton  solution  of  (3)  In  order  to 
define  the  Inviscid  velocitless  In  (11)  and  (12).  The  expression  for  the 
total  rate  of  energy  loss  per  width  of  tank  of  a  KdV  internal  soli  ton  Is 
complicated,  but  for  A  0, 


where  li  Is  the  maximum  wave  amplitude,  W  Is  the  tank  width  and  we  have 
assumed  vi  •  v2  ”  v  (because  a  «  1).  Of  the  three  terms  added  together  In 
(13),  the  first  represents  the  energy  loss  from  the  lower  fluid  due  to  the 
boundary  layers  at  the  tank  walls,  the  second  represents  a  similar  loss  from 
the  upper  fluid,  and  the  third  represents  the  energy  loss  at  the  Interface. 
In  the  particular  set  of  experiments  In  which  we  will  test  this  theory 
(hi  •  45  cm,  h2  *  5  cm,  W  ■  39.4  cm),  about  1/3  of  the  total  energy  loss 
occurs  at  the  Interface. 

The  total  energy  per  unit  width  of  tank  for  an  Internal  KdV  sollton 
with  A  «  1  Is  (from  Keulegan,  1953) 
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8p2  Aghjh^J^I^ 


(14) 


and  the  energy  loss  given  by  (13)  comes  entirely  from  this  source.  Between 
(13)  and  (14),  one  obtains  an  equation  for  the  slow  decay  of  a  KdV  soliton 
as  it  propagates  over  long  distances: 


\  hih2  J  y  hih2  j  ThT^IP^ 


(15a) 


Here  tio  is  the  initial  wave  amplitude,  and  n  is  its  amplitude  after  it  has 
travelled  a  distance  S.  If  a  <<  1> 
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12(gA)^'‘*  (hi  +  h2)^'‘‘ 
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W  2h2h2 


(15b) 


in  the  general  case  (0  <  a  <  1), 
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Lhj  +  h2 


1/2 


hi 


/  g/ J 


(15c) 


where  c^/g  is  the  smaller  root  of  (1)  and  (p.y)  are  given  below  (3).  In  the 
limit  A  0,  c^/g  ♦  (Ahih2)/(hi  +  h2),  and  (15c)  reduces  to  (15b). 
Moreover,  in  the  three  limits  given  in  (5),  (15)  reduces  to  the  results  of 
Keulegan  (1948)  for  the  viscous  decay  of  a  surface  solitary  wave. 

Koop  and  Butler  (1981)  used  an  approximate  version  of  (15)  that  amounts 
to  replacing  (7b)  with 


Uj  ♦  0  ,  U2  ♦  0  as  2  0  .  (16) 

For  a  <<  1,  their  method  gives  the  correct  answer  if  hi  «  h2,  but 
overestimates  the  damping  due  to  the  interfacial  boundary  layer  for  any 
hi  y  h2.  The  special  case  in  which  their  method  gives  the  correct  decay 
rate  (hi  *  h2)  happens  to  be  the  configuration  in  which  the  coefficient  of 
the  nonlinear  term  is  (3)  vanishes,  so  that  there  is  no  solitary  wave. 
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2  Comparison  with  Experiments 


The  experimental  data  that  we  will  use  to  test  this  theory  are  those 
shown  in  Figure  2  of  I.  Recall  that  in  those  experiments,  a  piston  at  one 
end  of  a  long  wave  tank  generated  . both  surface  and  internal  waves,  which 
then  propagated  down  the  tank.  A  vertical  plate  18.8  m  from  the  downstream 
edge  of  the  piston  was  carefully  lowered  into  the  water  after  the  surface 
waves  had  passed,  so  that  the  much  slower  internal  waves  could  be  measured 

without  disturbance.  The  fluid  depths  were  hi  «  45  cm,  h2  *  5  cm,  and  the 

density  difference  was  t  •  0.’^8. 

Based  on  Figure  3  of  I,  the  lead  waves  measured  at  each  of  the  last 
five  stations  may  be  considered  a  KdV  soliton,  to  within  experimental  error. 
Consequently,  we  may  use  the  lead  wave  at  x/(hih2)l/2  >  33.3  as  initial  data 
for  (15),  which  then  predicts  the  amplitude  of  the  lead  wave  at  the  next 
four  stations.  As  shown  in  Figure  la,  (15)  predicts  far  more  decay  than  is 
observed.  Whereas  French  (1969)  found  that  Keulegan's  (1948)  formula  for 
surface  waves  predicted  the  attenuation  to  within  about  12%,  errors  up  to 
30%  are  shown  in  Figure  la. 

Another  way  to  compare  (15)  with  these  data  is  shown  in  Figure  lb. 
Here  we  consider  the  lead  waves  measured  at  x/(hih2)l/2  «  33.3,  60,  100,  and 

151  to  be  initial  data  for  (15)  in  four  separate  tests,  and  ask  for  the  wave 

amplitude  predicted  at  x/(hih2)l/2  ■  191.  In  this  graph,  the  observed  wave 
amplitude  is  always  the  same  Cn  *  0.47  cm).  Here  (15)  predicts  too  much 
decay  for  long  propagation  distances  and  too  little  for  short  distances! 
Clearly  (15)  by  Itself  Is  Inadequate  to  explain  the  observed  attenuation  of 
Internal  solitary  waves,  even  though  It  Is  directly  analogous  to  the  much 
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more  successful  formula  of  Keulegan  for  surface  solitary  waves. 

We  propose  the  following  conjecture  to  resolve  this  paradox.  The  lead 
wave  in  Figure  2  of  I  was  not  propagating  into  a  quiescent  fluid,  as 
required  by  (15).  It  was  preceded  by  the  faster  surface  wave  (not  shown  in 
the  figure).  As  the  surface  wave  propagated,  it  lost  energy  both  by  viscous 
dissipation  and  by  leaving  kinetic  energy  in  the  boundary  layers  it 
generated.  This  residual  kinetic  energy  subsequently  was  lost  by  the 
boundary  layers  only  on  a  (very  slow)  diffusive  time  scale.  Therefore  the 
fluid  into  which  the  internal  wave  propagates  was  not  quiescent,  but  had  a 
mean  motion  down  the  tank.  Following  this  reasoning,  the  internal  soliton 
should  have  decayed  less  than  predicted  by  (15)  while  it  propagated  down  the 
tank  and  more  than  predicted  by  (15)  as  it  propagated  back  up  the  tank  (and 
met  not  only  the  residual  boundary  layer  of  the  faster  surface  wave  but  its 
own  residual  boundary  layers  as  well).  The  relative  position  of  the 
reflecting  wall  is  shown  explicitly  in  Figure  1,  and  all  of  the  data  are  in 
qualitative  agreement  with  this  hypothesis. 

Thus,  we  conjecture  that  (15),  which  requires  that  the  internal  soliton 
propagate  into  a  quiescent  fluid,  does  not  predict  the  observed  decay  of  the 
Internal  solitons  in  Figure  2  of  1  because  the  residual  boundary  layers  from 
the  faster  surface  wave  are  still  active.  Unfortunately,  this  conjecture 
leaves  us  without  any  means  to  predict  internal  soliton  amplitude 
realistically.  However,  it  does  emphasize  the  importance  of  these  residual 
boundary  layers  in  wave  tank  experiments  on  solitary  waves.  In  particular, 
attempts  to  measure  experimentally  the  phase  shifts  of  solitons  due  to  their 
Interaction  probably  should  be  regarded  as  inconclusive  until  these  residual 
boundary  layers  have  been  taken  into  account  properly. 
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Figure  Captions 


1.  Decay  of  amplitude  of  the  lead  sollton  of  Figure  2  of  I  as  a  function  of 
propagation  distances:  (a)  Distance  measured  from  Initial  location  at 

x/(hjh2)w2  ■  33.3:  measured;  _ _  predicted  by  (15); 

(b)  Distance  measured  from  final  location  at  •  191'* 

measured;  0.  predicted  by  (15). 
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